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Abstract. Dispersive averaging effects are used to show that KdV 
equation with periodic boundary conditions possesses high frequency 
solutions which behave nearly linearly. Numerical simulations are pre- 
sented which indicate high accuracy of this approximation. Further- 
more, this result is applied to shallow water wave dynamics in the limit 
of KdV approximation, which is obtained by asymptotic analysis in 
combination with numerical simulations of KdV. 



1. Introduction 

The study of the dynamics of high frequency waves, see e.g. [9], has 
been motivated by the so-called quasilinear phenomenon in optical commu- 
nication, where it was observed that spatially localized pulses evolve nearly 
linearly. The dynamics of high frequency waves can also be motivated by the 
well-posedness results in spaces of low regularity for dispersive PDEs, see 
e.g. [HE1I6]. These papers indicated that a subtle high frequency averaging 
effect took place in the nonlinear dispersive dynamics making these results 
possible. 

More recently, in [8| and [3], KdV was studied with regard to this aver- 
aging effect. In [8], near-linear dynamics was established for high frequency 
initial data and in [3] a new elegant proof of well-posedness in H s , s > was 
found using explicitly high frequency averaging effects. 

The purpose of this article is twofold. First we establish near-linear dy- 
namics in KdV under weaker and more natural assumptions than [8j. The 
proof relies on the so-called differentiation by parts technique (which is a 
variant of the normal form procedure) from [3] . Secondly, we investigate how 
near-linear dynamics for KdV can be extended to the water waves problem. 
We use the standard derivation of KdV in the long wave, shallow water 
approximation to obtain physical parameters for which near-linear behavior 
might be observed. 

We should note that for KdV on the torus or a circle the linear solution 
is periodic in space and time and thus one does not have dispersive decay. 
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It is also expected that the solutions of KdV on the torus will not be ap- 
proximated by the linear evolution for infinite time. Therefore the proof 
that the nonlinear evolution is almost linear (in the sense of the subsequent 
Theorem) on a finite but large time scale in a way provides evidence that dis- 
persion phenomena are not completely muted on the periodic setting. The 
seminal papers of Bourgain, [H E] by establishing Strichartz type estimates 
for periodic dispersive equations was probably the first step towards these 
new developments. 

Theorem 1. Consider the real valued zero mean solution of the KdV equa- 
tion 



Note that the difference between the actual solution and the solution of 
the linear KdV is small in L 2 for t up to e~ 1+ . 

We note that near-linear dynamics for high frequency solutions is easier 
to establish on unbounded domains, as the solution disperses to infinity 
and weakly nonlinear theories apply. On bounded domains (i.e. with 
periodic boundary conditions) the solutions cannot scatter to infinity. For 
the NLS in the 2d torus, Colliander et al, [7], gave recently a nice proof. 
Theorem [T] is an intermediate result between the iterated linear solutions 
and the situation at infinity. The nonlinearity averages out since dispersion 
will cause high harmonics to oscillate rapidly. Throughout this paper we 
assume that we have global well-posed solutions which have additional 
regularity properties. The interested reader should see [5] for the details. 
In addition, smooth solutions of KdV satisfy momentum conservation: 



Because of the momentum conservation law we can modify the equation 
adding a harmless term and thus only consider mean zero solution. This 
will imply the the Fourier series representation for the solution will have 
nonzero Fourier modes, an assumption that we consistently make in our 
paper. In particular notice that all the norms are restricted to this subclass 
of smooth solutions. In addition we use the conservation of energy, 



Uf — U xxx -\- UU X 

on T x IR with the initial data u(x, 0) = 4>(x) satisfying 

U\\2 = l, ||0||tf-l/2 =E«1. 

Then, for each t > and small 5 > 0, we have 






The KdV equation is locally well-posed (LWP) in L 2 (T), [5]. Due to energy 
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conservation KdV is globally well-posed and u G C(M; L 2 (T)). Kenig, Ponce, 
Vega, [12) . improved Bourgain's result and showed that the solution of the 
KdV is LWP in H S (T) for any s > — k- Later, Colliander, Keel, Staffilani, 
Takaoka, Tao, [6], showed that the KdV is globally well-posed in H S (T) for 
any s > — | thus adding a LWP result for the endpoint s = — \. Recently 
T. Kappeler and P. Topalov, [11] extended the latter result and prove that 
the KdV is globally well-posed in H S (T) for any s > — 1. 

The main idea of our proof runs as follows. First write KdV 

on the Fourier side, 

dtu k = y ^2 u kl u k2 - ik 3 u k , u fc (0) = 4>{k). 

k 1 +k 2 =k 

Then, using the identity 

(k t + k 2 f -k\-kl = 3(A;i + k 2 )k x k 2 , 



and the transformation 

v k {t) = Uk {ty kH 



the equation can be written in the form 

d t v k = ^ Yl ^ kklk2t ^v k2 . 

k±+k2=k 

The substitution 

u k {t) = e- uk \ k (t) 

eliminates the linear term ik 3 which has the highest growth at infinity and 
introduces oscillating exponentials into the nonlinear term. We have to 
show that v k stays almost constant for large times under our high frequency 
assumption. We cannot neglect the averaging effects of the exponent. 
Without the exponential factor the above system corresponds to Burger's 
equation which exhibits strongly nonlinear dynamics. Informally speaking, 
we mainly have three types of terms: 



a) Low frequency harmonics, with n\^n 2 small, give negligible contri- 
butions because of the high frequency assumption. 
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b) Intermediate terms are few in number as the Diophantine equa- 
tion rixnin = N has few solutions. 

c) High frequency harmonics, with n\,n2 large, are well- averaged by 
the exponent. 

We already mentioned that the method we use was inspired by [3]. 
Originally it was developed by Babin, Mahalov and Nicolaenko, [IJ in 
studying the global regularity of solutions of 3D problems in hydrodynamics 
(Navier-Stokes or Boussinesq system). In their framework the presence 
of high-frequency waves lead to destructive interference and weakened 
the nonlinearity through time averaging allowing one to prove global 
regularity. For the KdV the high Fourier modes of the linear term generates 
high-frequency oscillations which make the nonlinearity milder. There is an 
analogous phenomenon with the propagation of regularity to the Burger's 
equation with fast rotation 



For large |0| the nonlinearity weakens and the life-span of the solution is 
prolonged. So large oscillations is what separates the bad behavior of the 
classical Burger's equation and the good behavior of the Burger's equation 
with fast rotation. The same method has recently been applied by Kwon 
and Oh, [14] . to prove unconditional well-posedness for the modified KdV. 
The method of differentiation by parts helps to establish a priori estimates 
only in the C^H^. norms for any s > \. This is the heart of the matter 
in proving unconditional uniqueness, that is uniqueness of solutions to the 
modified KdV equation in the space C^H^ alone. 

The second motivation for our work comes from the various connections of 
the high frequency averaging process that we describe with certain aspects 
of the water wave theory. The dynamics of surface water waves has been 
an important object of study in science for over a century. Soliton solutions 
and integrability in P.D.E.'s are two examples of remarkable discoveries 
that were made by investigating water wave dynamics in shallow waters. In 
more recent times, the so-called rogue waves have been under an intense 
investigation, see for example [13} [T51 [2U] and the references therein. These 
unusually large waves have been observed in various parts of the ocean in 
both deep, see e.g. [16], and shallow water, see e.g. [IB] , motivating scientists 
to suggest various mechanisms for rogue wave formation. 



ut + uu x = iQu u(x, 0) = <p(x). 



Using Duhamel's formula the solution can be written as 
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In the case of shallow water, one normally does not work with the full 
water wave equation but uses approximate models to study the evolution, in 
particular the formation of rogue waves. These models are nonlinear disper- 
sive equations such as KdV, Boussinesq approximations, etc. In particular, 
KdV describes unidirectional small amplitude long waves on fluid surface. 
See, e.g. [T7] for applications of KdV to rouge waves in shallow water. 
Since rouge waves correspond to concentration of energy on small domains, 
one might argue that higher frequencies play important role in rouge waves 
formation. 

Here we provide some evidence, based on asymptotic expansions and nu- 
merical simulations that for sufficiently high frequency initial data, one- 
dimensional spatially periodic surface waves in shallow water exhibit near- 
linear behavior. Thus, linear theories of rogue wave formations can be ex- 
tended to nonlinear high frequency regime. 

Clearly, one has to be careful when considering short wave solutions for 
the equations obtained in the long wave approximations such as KdV. How- 
ever, we show that there is a set of parameters when our high frequency 
solutions correspond to a realistic physical scenario in shallow water waves, 
see Section [H 

2. Normal form reduction using "differentiation by parts" 

In this section we apply a variant of normal form reduction, called differ- 
entiation by parts [5J, to bring the equation to a more convenient form in 
which low order resonant terms are separate from the other terms. Using 
the Fourier series representation 

u (x,t) = Y,Mty kx 

with 

1 f n 

u k := u(k) = — / u(t,x)e~ lkx dx 
27 r Jo 

we express KdV as an infinite system of ordinary differential equations 
d t u k = y ^2 u ki u k 2 ~ ik 3 u k , u k (0) = 4>(k). 

k±+k2=k 

Notice that since the solution is real valued we have that u k = U- k - 
Changing the variable 

v k (t)=u k (t)e ikH 
(notice again that v k = v_ k ), and using the identity 



(i 
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(fci + k 2 f -k\-kl = 3(fci + A: 2 )A:ifc2, 

we obtain 

2 

Since e i3fcfc i fc 2t = e i3fcfcifc 2 t) differentiation by parts and ([I]) yields 



(1) ^fc = y ^ eW2 V^ 2 , « fc (0) = ^(fc). 

fcl+fc 2 =fc 



(1 v v e 3ikk 1 k 2 t Vk \ 

2^ E 3^2 2 ~2* E k^ 9tK ^ 
fcl+fc 2 =fc / fci+fc 2 — 



k 



1 / \ e 3ikkik 2 t Vk \ ^ ^ e 3ifcfcifc 2 t 

\fci+fc 2 =fc / fci+fc 2 =fc 

Note that since do = 0, in the sums above A;i and £?2 are not zero. The 
last two terms are symmetric with respect to k\ and k 2 and thus we can 
consider only one of them. Using Q we have 



g3ifcfcifc 2 t g3ikk±k2t I 

E -^-^^ = 2* E M E 

k\+k 2 =k k=k\+k 2 \u+A=fc 2 



e 3itfc 2M A 



1. Vk^nVX 3it]kk 1 (n+X)+n\( f J,+X)] 

2 ^ hi 

k=k±+fj,+\ 

We note that \i + A can not be zero since + \ = k 2 . Using the identity 

kk\ + fi\ = (ki + fi + A)&i + /xA = (fei + /u)(A;i + A) 
and thus by renaming the variables k 2 = /i, k% = A, we have that 

e 3ifcfcifc 2 i ^ e 3it(fci+fc 2 )(fc 2 +fc 3 )(A;3+fci) 

E h,h n Vk l 9tVk 2 = 9* E V kl V k2 V k3 . 

fc 2 +fc3^0 



hk 2 Kl 1 K2 2 ^ k 

fcl+fc 2 = fc k 1 +k 2 + k 3 =k 



All in all we have that 



<h [ v k - -B 2 (v,v) k j = ^R 3 (v,v,v) k 
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where 

B 2 (u,v) k = ^ ^ 

k 1 +k 2 =k 

and 

e 3it(k 1 +k2)(k2+k 3 )(k3+k 1 ) 



R 3 (u,v,w) k = 



fci 

kl+K2+K3=K 
k 2 +k 3 ^0 



-U kl V k2 Wk 3 . 



Now let's single out the terms (resonant terms) for which 
(2) ( kl + k2 )(k 3 + h) = 

and write 

R 3 (v, v, v) k = R 3r {v, v, v) k + R 3nr (v, v, v) k 

where the subscript r and nr stands for the resonant and non-resonant 
terms respectively. Thus 



R 3r (v,v,v) k = } y 



Ufcl«fc 2 «Jfe 3 



h 

k 2 +k 3 yt0 

and 

nr e 3it(k 1 +k 2 )(k 2 +k 3 )(k 3 +k 1 ) 



R 3nr (v,v,v) k = 



fel 

fci+fe2+A;3=fc 



-V kl Vk 2 Vk 3 . 



The set for which ([2]) holds is the disjoint union of the following 3 sets 
Si = {fei + k 2 = 0} n {k 3 + ki = 0} 4^ {fei = -fe, fe 2 = k, k 3 = k}, 

5 2 = {kt + k 2 = 0} n {fe 3 + fci / 0} ^ {fei = j, fe 2 = -j, fe 3 = k, \j\ + k}, 

5 3 = {fe 3 + fei = 0} n {fei + fe 2 / 0}} {fei = j, £; 2 = k, k 3 = -j, \j\ / k}. 



Thus 



s 
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R 3 r{v,v,v) k -2^2^ — r — - _ k + v >° Z> —j— +v k 2^ —— • 

\j\^k \j\^k 

Note that the second and third terms in the sum above are identically zero 
due to the symmetry relation j f-> —j. Thus 

t> I \ Vk \ |2 

We obtain 

dt (v k - ~B 2 (y,v) k J = -^v k \v k \ 2 + ~R 3nr (v,v,v) k . 

Since the exponent in the last term is not zero we can differentiate by parts 
one more time and obtain that 

nr e 3it(k 1 +k 2 )(k2+k a )(k a +k 1 ) 

R3nr(v,V,v) k = ^ 7 V kl V k2 V k3 = 

ki+k 2 +k 3 =k 1 

I I ™\ e 3it(k 1 +k2)(k2+k 3 )(k 3 +k 1 ) 

—d t B 3 (v,v,v) k - — V J-7-, — u — , ? w, , r v x 
3^ 3z fc + k2){k2 + k 3 )(k 3 + fci) 



{dtv kl v k2 v k3 + d t v k2 v kl v k3 + d t v ks v kl v k2 ) 

where 

™; e 3jt(fci+fc 2 )(fc2+fc 3 )(fc3+fci) 
S3(n ' U '" )fc = fe+ S k h(k 1 + k 2 )(k 2 + k 3 )(k 3 + kl ) Uk ^ Wk - 

As before we express time derivatives using Q. The terms containing dtv k2 
and 9tWfc 3 produce the same expressions and a calculation reveals that 



e 3ii(fci+fc 2 )(fc 2 +fc 3 )(fe3+fci) 

■X 



(d t v kl v k2 v k3 + d t v k , 2 v kl v k3 + d t v kz v kl v k2 ) = iB A {v, v, v, v) k 

where 

Bi(u, v, w, z) k = ^B\(u, v, w, z) k + Bf (it, v, w, z) k . 
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From now on J2* means that the sum is over all indices for which the 
denominator do not vanish. The term corresponding to dfV kl is 



Bl(u,v,w,z) k = 



e 3iti/)(fci,fc 2 ,fc3,/c 4 ) 



UkiVk2W k3 Z ki , 



ki+k2+k3+k4=k (*1 + **)(*! + k 3 + k 4 )(k 2 + k 3 + k 4 ) 

and the sum of the terms corresponding to dtv k2 and dtv k3 is 

* e 3itf(ki,k 2 ,k3M)(k 3 + k 4 ) 

B 4 (u,v,w,z) k = > — — — — — — —r-ru kl v k2 w k3 z ki . 

kl+k2 f^ +k4=k k i( k i + Mfa + h + k A )(k 2 + k 3 + fc 4 ) 

The phase function ^ will be irrelevant for our calculations since it is 
going to be estimated out by taking absolute values inside the sums. For 
completeness we note that it can be expressed as 

(fci + k 2 + k 3 + k A ) 2 -k\-kl-kl- fcf . 

Hence for R 3nr (v, v, v) k we have: 

R 3n r{v,v,v) k = ^-.d t B 3 (v,v,v) k - ^ (^B\{v ,v ,v ,v) k + Bj(v,v,v,v)k^j . 

If we put everything together and combining the two B4 terms in one we 
obtain 

(3) dt(v k - ^B 2 (v,v) k - ^B 3 (v,v,v) k ) = + ^B 4 (v,v,v,v) k , 



where 



^ (t,)fc= E k^ 2 

ki+k2=k 



* pj3(fcl+fc 2 )(fcl+fc3)(fc2+fc3)t,,, ,„ ,„ 
B (v) =- \ fc l fca fc 3 

fcl(fcl+fa)(*l+*s)(fe + *s) 



! * e #(fci,fc2,fc 3 ,fc4)t(2fc 3 + 2fe 4 + feiKw fc2 w fc 3W fc4 



2 , 1+ , 2 ^ +fe4=fc + fc2)(fcl + *S + fc4)(ib2 + ifes + ifc4) ' 

where if)(ki, k 2 , k 3 , k±) = (k\ + k 2 + k 3 + /C4) 3 — kf — k% — k\. 
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3. Proofs 

Notation: To avoid the use of multiple constants, we write A < B to denote 
that there is an absolute constant C such that A < CB. We will also use 
frequently the notation A < B(rj—) if for any 7 > 0, A < C^B(j] — 7). 
Similar notation will be used for A < B(rj+). Finally, for s £ 1R, we define 
the homogeneous Sobolev norm 

\w\\ H s = (Y J \k\ 2s w k \ 2 ) 1/2 . 

We have 

Proposition 1. The following a-priori estimates hold 



(4) \\B 2 (v)\\ H „ 1/2 <\\B 2 (v)\\ L 2 <\\v\\^_ 1/2 , 

(5) 11^3(^)11^-1/2 < \\B 3 (v)\\ L 2 < \\v\\%_ 1/a \\v\\v, 

(6) \\B4(v)\\ L 2< \\v\\^_ 1/2 \\v\\li, 

(7) \\B 4 (v)\\ A . 1/a < \\vfr_ 1/2 \\v\\lt, 

(8) Ib&AII^ ^ IMI#-i/aNljy"- 



Now we will prove Theorem [T] using Proposition [T] 
Proof of Theorem [II First note that 

\\u(-,t)-e td *(t)\\ 2 = \\u k (t) -e ifc3 *n fe (0)||^ (fc) = \\v k (t) - v k (0)\\ e2{k) . 

We will estimate the R.H.S. using ([3]). Integrating ([3]) from to T we have 

v k (T) - v k (0) = \b 2 {v(T)) - ~B 2 («(0)) + ^B 3 (v(T)) - ^B 3 (v(0)) 

6 J k 18 J 41 ,k 

The estimates in Proposition [TJ and the fact that for each t, \\v(t)\\ L 2 < 1, 
imply that 

(9) \\v(T) - v(0)\\ L 2 < \\v(T)\\%_ 1/2 + ||«(0)||^_ 1/a 

(10) 
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(11) \\v(T) - v(0)\\h-V> < \\v(T)\\%_ 1/2 + |K0)||^_ 1/2 
(12) 



+ 



''(')l|^_ 1/2 + |K./|l#-l/a 



)dt. 



Since ||v(0)||^_i/2 = ||(/>|| i j_i/2 < e, the inequality (11) and the continuity 
of the solution in L 2 (and hence in H^ 1 / 2 ) imply that, for all T < £~ 1+ , 
||v(r)||jj_i/2 < e. Using this in ^ implies that for T 



\\v(T)-v(0)\\ L2 <e 2 + Te 1 -. 

Proof of Proposition [J[ We start with ^ . Using vq = 

\\vl/k\\ e 2 < \\vk/y/k\\j<x>\\v k \\p < \\vk/y/k\\p\\v k \\p = \\v\\ 2 jj_ 1/2 \\v\\ L 2. 
We continue with Q. It suffices to estimate B 2 in L 2 : 
\\B 2 (v) 



□ 



Ml £ 

ki+k 2 =k 



k\ k 



\K2 



< II N N 11 
- II fc fc \\p 



\wm\%/3 < 1^11^-1/2, 



where we used the Young and Holder inequalities in the second and third 
inequalities respectively. Now consider Q: 



E 



k\+k 2 +kz=k 



e i3{k 1 +k 2 )(k 1 +k 3 )(k 2 +k s )t v ^ v ^ Vk 2 



h{h + k 2 )(h + k 3 )(k 2 + k 3 ) 



< 



E 

^1+^2+^3=^ 



W2 



M \v k2 



+ k 2 \\h + k 3 \\k 2 + k 3 \ ^\kl\ y/\k 2 
By Cauchy Schwarz we estimate this by 

N 



kfca 



<2 



E E 



E 



l^nil I^ri2l I 



J n 3 



< / N 

" SUP V ^ |A;i||A;i + fc2| 2 |A;i + fc3| 2 |fc2 + ^l 2 



fcl+fc2+fc 3 =fc 



* I |2 I |2 

E \ v ni\ \ v n 2 \ 1 12 



■■•1-/ •- ' — * Inil \n 2 

ni,n 2 ,n 3 1 11 1 z 
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II. .1,2 V- N 



|4 II ||2 

\h-^\\ v Wv su p 2^ 



It remains to show that the supremum above is finite. Note that the supre- 

< 



sup ]T 



N < \ h 



where we used the fact that, for ki,k 2 / k, \k — ki\\k — k 2 \ > \k\ — k 2 \- 
Now to estimate this sum, consider the cases \k\\ > 2\k 2 \, \k\\ < \k 2 \/2, and 
I /ci| ps \k2\ separately. In the first case, the sum is < 



|fcl|>2|fc 2 | 

In the second case, we have 



\k 2 \ 2 \h 



< 00. 



V 1 < ^ iogCND 

^ \k 2 \ z \k x \~ ^ \k 2 f 

|fei|<|fe 2 |/2 1 1 1 1 1 1 



In the third case we have 

* 1 * 1 

Finally, we consider B4. First note that 

in ( v \ 1 < v ^wvy 

' 41 ~ 1*1 + fe||fcl + *S + fc4||fe + fcs + fc4| 



^1+^2+^3+^4=^ 



\- h^Wg^fcjl _. R 1 W , d2/ \ 
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First we consider the L 2 norm of B\. Applying Cauchy Schwarz as in the 
case of B3 , we have 

\\Bl(v)\\bZ 

* 1 

^ \ki\*~\vk 4 \ Ki^fa^fcgJ 

kl+k2+k3+ki=k \h + k 2 \\h + k 3 + k^-\k 2 + h + k 4 \ \h\2-\k1 + k 3 + /c 4 |2" 

, A l^il 1 "!^! 2 

< sup 2^ 



*\+ fc2 ^=* + fc2 l 2 l fcl + * 3 + fc *l 1 -l fe + fc3 + fc4 ' 2 



* 1 1 2 

^ ^ I ''in ' ' 



r^i 1 1 Ini + 713 + n4| 1+ 

rai,n2,ri3,n4 111 1 



Note that the sum in the parenthesis is < || v||^._ 1/2 _ l _ |M|^ 2 summing in 
714 first). We estimate the supremum by eliminating k 3 in the sum as follows 

m \ l^il 1 '!^! 2 

I, ||2 l^ll 1 ' 

- \\v\\ L 2 sup ^ + _ _ 

Using \ki + k 2 \ 2 \k - k 2 \ 1 ~\k - ki\ 2 > |fci| 1_ |A;i + k 2 \ 1+ \k - ki\ l+ we have 



The last inequality follows by summing first in k 2 then in k\. Now consider 
the L 2 norm of B\. Similarly, we obtain 

|| R 2/^||2 <am bfci| 2 |^fc 2 | 2 kfc 3 | 2 

\\B 4 (v)\\ L 2 < sup ^ T^n= 

^1+^2+^3+^4=*: 



y I 2 

I 



|ni| 1+ |rii + n 2 \ 2 \n 2 + n 3 + n 4 | 2 / ' 

The sum in parenthesis is < 11^11^2 by summing first in 723, then n^, then 
n 2 , and then in n\. Finally 

V- l^i| 2 |^ 2 | 2 bfc 3 | 2 . V- bfci| 2 |^ 2 | 2 bfc 3 | 2 _ ,|4 I. ,||2 

sup 2^ ^ r^p - IMMMI j-w 

^1+^2+^3+^4=*: ki,k 2 ,k 3 

Combining the estimates for B\ and B 2 , we obtain ([6]): 

II^4^)||l2 < HwlliallwII^-i/JH < IMlitlMI^-i/2- 
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It remains to prove ([7]). We start by estimating B\. Applying Cauchy 
Schwarz as above we have 

ii Di/ \m2 < Ifcil 1 "!^! 1 "!^! 2 

wm^n - r ki+i ± +i4 = t \m + wih + h + h nt 2 + k3 + w 

* I 1 2 

| ^ni ^Tl2 ^Tl3 I \ 

^ Ini 1 1— llT-2 1 1_ 1^1 + ?T-3 + 77, 4 | 1+ / ' 

Note that the sum in parenthesis < IMlia IMI^-_i/2+- Eliminating ks in the 
first sum we have 



\Bl(y)\\ 2 6 _ 1/2 < |Mli 2 |M|^_ 1/2+ sup 



Ifcil 1 -!^! 1 -!^! 2 



k h Tu \k\\k 1 + k 2 \ 2 \k-h 2 \ 1 -\k-k 1 \ 



sr Ifcil 1 "!^! 1 - 

L 2\\o U _ 1/2+ sup ^ + _ _ 



II Il4 11 ||4 
< \\V \\t3 U" 



Note that the sum above is bounded by 



v i**r- 1 y m i - 

2^1 \h, 4- fcJ2|fc _ fcoH-lfc _ h, |2 + 2-> 



, , Ifci + ^PIfc-fcali-lfc-A;!! 2 ' |A;p 1 + A: 2 | 2 |fe-A :2 | 1 -|A:-A; 1 | 1 + 



- ^\k+k 2 \ 2 \k-k 2 \ i - + ^ \k\\k + fc 2 |i+ifc - fc 2 |i- ~ 4^ < 00 

In the first inequality we used (for a > b > 1) 



1 



< 1 



n + m| a |n| 6 ~ \m\ b ' 



which follows by considering the cases \n\ < \m\/2 and \n\ > \m\/2 sep- 
arately. In the second inequality we used \k + k 2 \\k — k 2 \ > \k 2 \ and 
\k\\k - k 2 \ > \k 2 \. Similarly, 



\Bj(v)\\l_ 1/2 <sup Yl 



|Wfcl| 2 |v*2| 2 |Wfc 3 



fc ^ l^il 1 -!^! 1 - 

ft fci+fc 2 +fc3+fc4=A: 11 11 

* * 111 I I o 

\ - \ - P 2 | K 4 | 

^ , \k\\ni\ 1+ \ni + n 2 \ 2 \n 2 + n% + nA 2 

k ni+n2+rt3+?i4=fe 
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Note that the supremum is < |MIl2|MI^--i/2+- Eliminating n.3 in the paren- 
thesis we obtain 



* 1 1 1— 1 1 2 

i R 2/ x 11 2 <■ 11 i|2 || ||4 l ra 2| Pn 4 l 

I^WIIff-i/a % \\v\\ L2 \\v\\^_ 1/2+ 2. | fc || ni |l + | ni+n2 |2| fc _ ni |2 

1- 



< 



I Il4 11 n4 n 2 



-1/2+ XT 



£2|| U || H _ 1/a+ ^ Iife||ni| 1 +|ni+n2| 2 |fc-ni| 2 

ni,n2,A: 



Applying (14) to the sum in k, we have 



* 1 1 1— 

in2/ mi 2 ^ 11 n4 11 11 4 r^2 ^11 11 4 11 n4 

|S 4 (v)||h-i/2 £ \H L 2\\v\\ H _ 1/2+ ^ I 12+1 — 12 % \Hi?\Hh-W 

ni,n 2 1 1 1 1 



The last inequality follows by applying (14) to the sum in ri\ and then 
summing in 712- This yields ((7]). □ 



4. Numerical simulations demonstrating high accuracy of 

approximation 

In this section, we present numerical evidence that for the initial data with 
sufficiently high frequency, linear KdV approximates very well nonlinear 
KdV. As the initial condition we use the first Hermit function, appropriately 
scaled 

The Figures [I] and [2] show the initial and evolved waves in KdV with 
periodic boundary conditions u(x + 2n) = u(x) for the time interval T = 2tt. 
Here e = 0.1 and the time step is At = 10 -7 . Note that the linear KdV 
evolution is 27r-periodic in time. Therefore, if the near-linear dynamics takes 
place, we should see nearly perfect return of the evolved data to the original 
profile. Both figures confirm such behavior. 

The Figure [3] demonstrates an obvious but important property that away 
from t = 2ttN, the evolved data is very far from the initial data. 
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Figure 1. Spectral data for initial and evolved waves. 
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Figure 2. The initial and evolved data. The time interval 
is 2tt, which is a period of linear KdV, therefore the initial 
and the evolved data are very close. 

5. Near-linear dynamics in water waves: asymptotics, 
numerical simulations and physical interpretation 

Rogue waves are large-amplitude waves appearing on the sea surface seem- 
ingly "from nowhere". Such abnormal waves have been also observed in 
shallow water and KdV has been used to explain this phenomenon [T7j. In 
the oceanographic literature, the following amplitude criterion for the rogue 
wave formation is generally used: the height of the rogue wave should exceed 
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FIGURE 3. Evolved data after short time. There is large 
distortion in physical space because of the strong dispersion. 

the significant wave height by a factor of 2-2.2 [16 1 . (Significant wave height 
is the average wave height of the one-third largest waves.) 
Major scenarios and explanations of rogue waves involve 

• probabilistic approach: rogue waves are considered as rare events in 
the framework of Rayleigh statistics 

• linear mechanism: dispersion enhancement (spatio-temporal focus- 
ing) 

• nonlinear mechanisms: in approximate models (e.g. KdV), for some 
special initial data, large amplitude waves can be created. 

Linear mechanism is very attractive as there are simple solutions leading 
to large amplitudes, while nonlinear mechanism requires rather special initial 
data, e.g. leading to the soliton formation. On the other hand, linear 
equations arise in the small amplitude limit which is too restrictive. 

Using near-linear dynamics in KdV one can experiment with another 
mechanism of large wave formation that combines linear and nonlinear de- 
terministic mechanisms. Our results indicate that for a special but rela- 
tively large set of initial data (characterized by the energy contained mostly 
in high frequency Fourier modes), the solutions of KdV equation behave 
near-linearly. It is then possible to construct large amplitude solutions using 
linear mechanisms of large wave formation. 

Here, we illustrate our approach with periodic boundary conditions. This 
model is not the most realistic one but appropriate to illustrate the concept. 

The KdV equation has been used to describe surface water waves in the 
small amplitude limit of long waves in shallow water. More precisely, two 
parameters are assumed to be small and equal 

amplitude / depth 
depth \ wavelength 

Our numerical simulations of KdV show that the near-linear dynamics phe- 
nomenon occurs when small parameter e characterizing high frequency limit 
(see the formula below), is only moderately small e = 0.4. The following 
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three figures show that with e = 0.4, there is still a clear presence of near- 
linear evolution of KdV for some reasonable time interval T = 1 . 

On the other hand, we will show that this value of e = 0.4 is sufficiently 
large so that KdV still approximates shallow water waves dynamics. 

As the initial data, we take the scaled 1st Hermit function 



u(x) 



4.5 



so that the energy J u 2 dx does not depend on e and is very close to 1. For 
numerical simulations, we use KdV in the form 

3 1 

Ut — ~UUx ~^U X xx 
I O 

as it appears in the derivation of KdV from the water wave equations (see 
below). Specific numerical parameters are: the length of periodic domain 
L = 2n. The number of modes M = 2 9 . Time step size At = 10" 7 with 
the time of the evolution T = 1. The discretization in space is given by h = 
L/M. We used the so-called Fornberg-Whitham scheme which is described 
in [10]. 

Now, using standard derivation of KdV from water waves equations, we 
recall the relation between physical parameters and rescaled dimensionless 
variables, see [TJ5], Chapter 13.11. 

Let ho be the depth when the water is at rest and let Y = ho + r] be the 
free surface of the water. Let a be a characteristic amplitude and / be a 





100 150 200 250 



Figure 4. Initial and evolved waves in KdV in Fourier space 
(right) and physical space (left). The nonlinearly evolved 
data in physical space is pulled back with reverse linear evo- 
lution, e~ Lt u{x,t), for proper comparison. Abscissa shows 
the number of Fourier harmonic (right) and spatial coordi- 
nate (left), while ordinate is the amplitude. The time of 
evolution is T = 1, e = 0.4. 
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Figure 5. The initial data and the solution after t = 0.2. 
While, e = 0.4 is not so small, the dispersion is sufficiently 
strong so after a short time the initial wave disperses over 
the whole periodic domain. 



characteristic wave length. Assume that 

a n hr. 

h I 2 

Both a and f3 are small parameters in the problem and they must be of the 
same order. 

Next, use the following natural normalization 

x = Ix, Y = h$Y, t = U/cq, Tj = arj, 

where primed variables are the original ones and Co = \fgh~Q. 

The formal asymptotic expansion leads to KdV with higher order correc- 
tions 

3 1 

Vt + rix + ^arjrjx + -(3rj xxx + 0(a 2 + (3 2 ) = 0. 

2 D 

Let X = x — t and T = at, so the equation becomes 

(15) Tfr + -rmx + Ivxxx + 0(a + P 2 /a) = 0. 

2 D 

One should expect that this approximation has accuracy of the order 0(a) 
for finite time T = 0(1), which implies t ~ a -1 and t' ~ 1/{cqq). 

Finally, since we modify our solution with another parameter e, we ver- 
ify that KdV approximation will still make sense for some choice of the 
parameters. 

First, let a = ft = 5 <C 1. Let us modify a and / with a e = -^a and l e = el 
which is consistent with our scaling of initial data. Then, we have 

_ Cle _ $ fl-^I-A 
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These are small with 5 = 0.01 and e = 0.4. On the other hand the "mis- 
match" in the equation (15) is 



a £ + 



+ 



-3.5 



Therefore, our high frequency regime may approximate water waves dynam- 
ics for example with the following parameters: a = 1 m, ho = 100 m, and 
I = 1000 m. 
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